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Abstract 

I consider an extension of General Relativity by an auxiliary non-dynamical 
dimension that enables our space-time to acquire an extrinsic curvature. Ob- 
tained gravitational equations, without or with a cosmological constant, have 
a selfaccelerated solution that is independent of the value of the cosmological 
constant, and can describe the cosmic speedup of the Universe as a geometric 
effect. Background evolution of the selfaccelerated solution is identical to that 
of ordinary de Sitter space. I show that linear perturbations on this solution 
describe either a massless graviton, or a massive graviton and a scalar, which 
are free of ghosts and tachyons for certain choices of boundary conditions. The 
obtained linearized expressions suggest that nonlinear interactions should, for 
certain boundary conditions, be strongly coupled, although this issue is not 
studied here. The full nonlinear Hamiltonian of the theory is shown to be 
positive for the selfaccelerated solution, while in general, it reduces to surface 
terms in our and auxiliary dimensions. 



1. Extension of General Relativity 

One simple way to parametrize the cosmic acceleration fT] is to introduce in the 
Lagrangian of General relativity (GR) the cosmological constant A ~ (10~^^ eV^. 
This is not quite satisfactory however, since the parameter A receives contributions 
from various scales of particle physics each of which is many orders of magnitude 
greater than (10^'^^ eV)"^. Without an underlying principle, cancellation between 
these contributions down to (10~^^ eV^ seems conceptually unlikely and technically 
unnatural [2j. 

Here we consider an extension of GR, such that for an arbitrary but given value of 
the parameter A, there exists a solution - requiring adjustment of certain boundary 
terms - that is independent of A. Furthermore, the observed cosmic acceleration 
will be due to a new parameter m with the dimensionality of mass, appearing in 
the extended GR Lagrangian. This parameter does not receive contributions from 
the particle physics; its value can be set to m ~ 10~^^ eV. The present approach 
does not explain the smallness of m; instead it gives a technically natural way of 
describing cosmic acceleration, with potential observational predictions that differ 
form those of GR with the cosmological constant. The present approach does not 
solve the cosmological constant problem either, but instead it reduces the problem 
to the choice of the boundary conditions in the classical gravitational equations, 
with everything else being quantized (more on this in section 2.) 

The gravitational field will be described by an extended metric tensor g^^{x,u), 
with /i, 1/ = 0, 1, 2, 3, which is labeled by a continuous dimensionless parameter u. 
The extended metric varies as g'^j^{x',u) = u!'^{x)uj^{x)ga/3{x,u), under the general 
coordinate transformations x'^ = [u!~^{x)]'^x'' . This leaves the extended interval 
dsf^ = gfj,i,{x,u)dx^dx'^ invariant. However, the matter fields do not depend on u; 
they will only couple to the metric tensor 

9tiA^) ^ 9f^uix,u = 0) , (1) 

with the relevant invariant interval being ds"^ = gfj,p{x)dx^dx'^ . 

Consider the Lagrangian density for the gravitational field (we use the conven- 
tions of [3] and also set Mpj = (SttGat)"^ = 2, unless stated otherwise) 

jC = ^R±m'j ^ du^g{k%-k'') , (2) 

where R is the Ricci scalar of the metric g^y{x), while k^y = ^dug^u, k = g^^k^u', 
all indexes in the Einstein-Hilbert term in ([2]) are raised by g^" , while those in 
the second term in ([2]) by g^^ . The value of k^y measures an extrinsic curvature 
of a (3 + l)-dimensional constant-w surface in certain coordinates in the "x — u 
space-time". The Lagrangian density ([2]) is covariant in (3 + l)-dimensions. 

We impose the Z2 symmetry on the fields g^v{x,u) = g^u{x,—u) across the 
hypersurface u = 0. Then, it is enough to consider the interval [0, 1] for the variable 
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u. Note that the "-u-dimension" is not dynamical since fields have no ordinary 
derivative terms there. Moreover, integration boundaries in u may in general take 
any finite value, which can be reduced back to the interval [—1, 1] by an appropriate 
rescaling of u and the parameter m, before specifying boundary conditions that 
could be sensitive to such rescaling. 

We refer to the surfaces u = 0, ±1 as fixed boundaries. Eq. ([1]) imposes one 
boundary condition on the M-dependence of the extended metric. This is not enough 
to determine completely the u-dependence of g^^{x,u), the second boundary condi- 
tion should also be specified. For this one can either impose the Neumann-type or 
Dirichlet-type condition at the either boundaries. For now we keep this condition 
unspecified and find various solutions that correspond to different choices of the 
second boundary condition. 



2. Equations, Solutions and Boundary Terms 

Let us start with the action of gravity plus "everything else": 

S = j d^x^{R + C{m,g)}±m^ j d^x ^ du>/^ {kl^ - P) . (3) 

Here (7) is the non-gravitational Lagrangian the fields in which couple univer- 
sally to the metric tensor g^i,{x), hence preserving the equivalence principle. We'll be 
looking at very low-energy phenomena (as compared to Planck's scale) and thus re- 
gard effective classical field describing large distance gravitational physics; 
thus, the gravitational part of the action will not be quantized (it can be regarded 
as the IPI effective action in which all the quantum loop effects are encoded in the 
coefficients of various terms). All the other interactions encoded in C{^,g) will be 
quantized. 

The Lagrangian ,g) will contain in general the cosmological constant gener- 
ated by particle physics. As noted earlier, it receives contributions from the scales of 
electromagnetic, strong and weak interactions; we denoted it hy K fund- Furthermore, 
quantum fluctuations of the non-gravitational flelds in C{^^g) will generate higher 
dimensional gravitational operators, such as B? , etc., which all are functions of 
g^y and are suppressed by the Planck's scale. Importantly, non of these terms, that 
are signiflcant in the UV, can change the effects of the second term in ([3D which 
switches on in the IR. Moreover, the second term in ([3]) does not get renormalized 
by the quantum loops of particle physics, since the particles couple only to g^u{,x) 
and cannot give rise to operators made of g^j_p. This can also be seen from the 5D 
representation of the model given in section 5. There, the matter flelds localized 
on the brane cannot renormalize the bulk terms because of geometric separation in 
extra dimension; the bulk terms stay unchanged, as long as gravity is considered to 
be a classical fleld theory with the effective IPI action. 
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The equations of motion obtained by varying the action ([3]) SS with the fixed 
boundary conditions in the u — space, Sg^^{x, u)\u=±i = (amended by the condition 
Sg^u{x,u)\x=boundary,u=o = wheu the bouudary is present in the x-space, in which 
case the Gibbons-Hawking (GH) boundary term |1] should also be introduced in the 
action) gives the following two equations for u = O"*" and < m < 1 respectively: 

G^, ± 2m' {k^, - g^,k) = T^J2 , (4) 

and 

^ {kg'-'' - kn] = h'^V^ {k' - kip) + 2v^ {k^% - k^^k) . (5) 



du 

Note that the right hand side (rhs) of Eq. ([5]) is traceless. 

Furthermore, equation (j4]) combined with the Bianchi identities implies that: 

D^'k^, = D,k , (6) 

where denotes the covariant derivative compatible with the metric. Eq. ([6]) 
should automatically be satisfied by any solution of (jlj). Note that Eqs. (jl]) and 
are similar to those of the DGP model [5] written in the 5D ADM [6J form (see, e.g., 
[7j). However, there are two significant differences: (a) What is the {55} equation 
in DGP is absent here; (b) In Eq. ([5]) there are no derivatives w.r.t. space-time 
coordinates, and thus it significantly differs from its DGP counterpart (what is the 
bulk {/ii^} equation). 

Equation ([5]) determines the evolution of the metric g^^ in the w-direction. This 
is a second order equation. One boundary conditions for it is specified by ([T]); 
pending the second boundary condition we find different dependence of the metric 
on u. The latter sets the value of the extrinsic curvature at m = 0"*", which by its 
turn determines 4D geometry via Eq. (jl]). 

We turn now to concrete solutions. In the absence of any matter stress-tensor or 
cosmological constant (T^,y = 0) the above system of equations has the Minkowski 
solution g^^{x,u) = t]^^ = diag{ — 1, 1, 1, l}^u- In general, for the class of extended 
metrics which are independent of u, the theory at hand reduces to GR. This would 
correspond to the choice of the boundary condition dug^u\u=o = 0, in addition to 
([1]). Thus, for g^u{x, u) = rj^j^iy + /i^,y(x) the fluctuations of the extended metric above 
the Minkowski solution describe a massless graviton. 

There exists a choice of the boundary conditions for which the linearized fluc- 
tuations describe a Minkowski space massive graviton; for instance, the Lagrangian 
([2]) with the minus sign in front of the second term describes the Pauli-Fierz massive 
graviton of (mass)^ = 2m^, with ^^,^(a:, m) = rj^u -l- (1 — \u\)h^i,{x) being a linearized 
solution selected by imposing the second boundary condition in the Dirichlet form: 

Hence, the theory ([2]) endowed with the appropriate boundary conditions gives 
a nonlinear completion of massive gravity. Remarkably, the Hamiltonian of this 
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theory does not suffer from tlie problem found in Ref. |8j in 4D massive gravity, as 
it will be shown in Section 4. Since Minkowski space is not a subject of a primary 
interest here, we will not elaborate on this branch of solutions further. 

Consider now a factorized expression for the extended metric 

g^uix,u) = a{u)gf,^{x) . (7) 

The rhs of Eq. ([5]) is identically zero for ([7]), and Eq. ([5]) reduces to d^a = 0. Hence, 
for M > we have a{u) = cq + ciu, where cq, ci are integration constants. The 
boundary condition ([1]), and ([7]) define the value of Cq = 1, while Ci has to be fixed 
by the second boundary condition. Below we consider various solutions that differ 
from each other by the choice of that boundary condition. 

For the second boundary condition specified in the following form 

dugf,u\u=o+ = TOi^uix) , (8) 

it is straightforward and not tedious to check that the system of equations (14]),([5]) 
admits a selfaccelerated solution: 

r^.ix, u) = {lT \u\)%u{x) , R{l) = 12m2 . (9) 

Here, 7^j/(a;) denotes the 4D de Sitter metric with the Hubble parameter H equal 
to m. This solution can describe the cosmic acceleration of the Universe, with the 
acceleration being due entirely to a geometric effect. In that regard, the growing 
solution in is similar to the selfaccelerated solution [21 [ID| on the DGP model, 
while the decaying solution to that of Refs. [11] . 

For the decaying solution in Qj the extended metric (jf^u vanishes at the bound- 
aries u = ±1, while the inverse of g^i, is singular, giving rise to a singularity of the 
extended Ricci tensor R made of g^^. However, since the "w-dimension" is nondy- 
namical, and all the matter and their interactions are located at m = 0, the extended 
Ricci tensor R evaluated at u = ±1 should not have a particular significance. More- 
over, this singularity is easily avoidable by changing in ([2]) the integration interval 
for u from [0,1] to [0,6], where 6 < 1 is some positive number. This would not 
change the equations (jl]) and ([5]) and the solution ([9]), but for any b ^ 1 one would 
need to add to the Lagrangian a surface term. The latter would guarantee that 
the effective Lagrangian obtained by integrating over the w-direction (i.e., by first 
substituting the metric ([7]) into the action and then varying it w.r.t. the metric 
g) gives the result consistent with the solution ([9]) obtained from the equations of 
motion (jl]) and ([5]). 

The Lagrangian with the surface terms for general b, which gives the selfaccel- 
erated solutions Qj, reads as follows: 

^ du^g {k% - + ^m^ + ^/^g\u=-b) , (10) 
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where C^"^^ = -3/(1^6). 

For the growing solution with the positive sign in iQ the above singularity is 
absent, however, even when 6 = 1 is chosen, for this solution one has to add to 
the Lagrangian the surface term given in flTU]) in order for the effective Lagrangian 
(obtained by integrating out the u-direction) to give the result consistent with the 
solution (Q that was derived from the equations of motion (jl]) and ([5]). Moreover, 
for the growing solution in (Q the surface term will be crucial for calculation of its 
energy in Section 4. 

Although we have constructed these surface terms via "inverse engineering" star- 
ing with the desired solutions, the straightforward statement is the following one: 
for the given boundary conditions and specified surface terms there are unique self- 
accelerated solutions corresponding to the two sign choices in ([2]). 

One could of course modify the second boundary condition ([8]) in various ways 
and obtain different solutions, to some of which we're turning now. 

For a nonzero homogeneous and isotropic stress-tensor there exists a solution for 
which the extended metric reads g^u{x,u) = (1 =i= i.^) ^ and the modified 

Friedmann equation in the standard notations takes the form 

H'-Cm' + - = ^-^p, (11) 

CI o 

where k = ±1,0 labels the 3D spatial curvature, and ( is an arbitrary integration 
constant that could be fixed only after imposing the boundary condition for e.g., 
gf,,.ix,u)\u=±i, or for 9„^^j.(x, m)|„=o,±i- 

If the stress-tensor contains the cosmological constant {SttGnT^i, = Afundd/iu) 
the value of ( can be chosen to cancel its contribution down to zero. This can be 
combined with the selfaccelerated solution obtained above. For instance, consider 
the Lagrangian with the cosmological constant and the choice of the plus sign in 
front of the second term in ([2]): 

Cb = ^{R - 2Ayw) +ni^ j du^ {kl, - + C^m^ (v^U=fe + v^U=-b) -(12) 

The corresponding equations (jl]) and ([5]) have a consistent solution: 

r^^ix, u) = il + C\u\ - \u\)%,{x) , R{^) = 12m' , (13) 

if = 3(C - 1)/C, where C = Afund/Sm' > 1. 

The result of this discussion is the following: for an arbitrary value of the cosmo- 
logical constant generated by particle physics A fund, one can choose the boundary 
conditions and surface term in (fT2|) such that the background solution describes 
an accelerated universe with the Hubble parameter that is independent of A fund, 
but instead is defined by the UV insensitive new mass scale m, introduced in the 
Lagrangian ([2]), or f|T2l) . 
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This scheme does not provide a dynamical mechanism for solving the cosmo- 
logical constant problem, as one has to adjust the boundary terms and conditions 
appropriately to get rid of Is. fund- However, it has an advantage over GR in the fol- 
lowing respect: GR, as well as the present model, at the classical level can entirely 
be specified by their equations of motion, without any reference to the action. The 
GR equations with the cosmo logical constant have no other solutions but the (A)dS 
solutions with curvature set by the value of the cosmological constant. In contrast 
with this, the equations of motion of the present theory with the cosmological con- 
stant do have solutions with curvatures that are not related to the cosmological 
constant. The above properties of the equations make no reference to the boundary 
terms. The latter come into the play only when the action functional is invoked. 

Hence, as long as gravity is treated classically while all the other interactions are 
quantized, the present approach reduces the cosmological constant problem to the 
choice of the boundary conditions in the classical gravitational equations. 

The fact that A fund can be removed by means of the boundary conditions which 
specify the otherwise arbitrary integration constant, is somewhat similar to what 
happens in the unimodular gravity [TB] where the cosmological constant can be 
fixed by superselection rules. However, a distinction between the two approaches 
is that the perturbations in the present case can be different from those of the 
unimodular gravity which are identical to the GR perturbations. 

In the context of inflationary cosmology, the present method would remove a 
constant piece from the inflationary potential, while retaining all the positive aspects 
of the slow-roll inflationary paradigm (note a similarity in this with Ref. |14j). 

As mentioned before, the theory ([2]) contains all the solutions of GR: using 
the factorized form ([7]) with a = 1 one would obtain just Einstein's equations for 
g^u- For the selfaccelerated universe a = 1 ^ \u\, and equation (jlj) for gf^,^ is the 
ordinary Einstein equation with the cosmological constant equal to 3m^. Thus, for 
instance, the dS-Schwarzschild solution of GR is also a factorized solution on the 
selfaccelerated background. Similar arguments apply to any other solution of the 
Einstein equations. Furthermore, there may well exist other solutions, e.g., for a 
static source, that do not have the factorized form ([7j). The latter would be selected 
from the factorized solutions by the boundary conditions. 

Factorized or not, the spectrum of linear and/or nonlinear perturbations about 
these solutions are determined by Eqs. (jl]), ([H]), which themselves may or may not 
have a factorized form ([7]), or depending on boundary conditions, could differ from 
the spectrum of GR. One example of this is the spectrum of linear perturbations on 
the selfaccelerated solution to which we turn in the next section. 
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3. Perturbations of the selfaccelerated solution 

We denote the deviation from the background metric as follows: 

g^lv{.x, u) = g^Jl^ix, u) + 5gf,^{x, u) . (14) 
Where, g'^'' is defined in (jH]). It is straightforward to derive that 



„1_ _ ^ _ - 2 

k^u = + Sk^j^ , k^y = ^—"^^u 1 k = k -\- Sk , k = g^ k^u = , (15) 



a 



where 



Sk,, = ^djg,, , 6k = j-^^^'^djg,, ± ^r^Sg,, ■ (16) 



An expansion of Eq. (jl]) in the linear approximation reads: 

SG,, ± 2m^ (5k,y - Sg^^k - g^Jk) = T,,/2 , (17) 

here T,, on the r.h.s. is the stress-tensor of a probe source which has nothing to do 
with the background; the variation of the Einstein tensor on the dS space is 

= -i {n6g,, - V^V'Sgau - V.V"^^?.^ + V^V.^^^) 

-l%u (V"V^5^7a/3 - nSg:) - 2HHg,, + \H^i,u6g^^ , (18) 

where V denotes the covariant derivative w.r.t. 7. The constraint ([6]), which in the 
linearized approximations reads 

± V^Sg,, + Wdjg,, = ±VJgl + V,9,<5^2 , (19) 

can be satisfied by the following gauge fixing condition 

V''5go.p = Vp5gZ. (20) 

Using the latter in equation f|T7|) . where we also substitute w? = H^, we obtain: 



- ^ (□%, - V,V,6g:) + 2HHg,, - \i,.HHg'^^ 

±H^ {djg,, - i^Mgl) = T,,I2 . (21) 

Taking trace of the above equation gives: 

T'iii'dJgl = Tl2. (22) 

One needs to solve equation ^ to obtain the u-dependence of the perturbations. 
For this one considers variation of its left and right hand sides separately at m > 0: 



^du^ip^'^T^dJgo.p T ^l^^'T^Sg^p - ^r"7^^5«5^?a/3 ± W'^l^'^Sg^^ . (23) 



Notice that all the equations presented above with the upper sign are equivalent to 
those with the lower sign provided that in the latter the replacement u — >■ —u is 
made. This will be reflected in our final solutions. 
The variation of the rhs of ([5]) equals to: 

l^^r^r^Sg^f^ ± Ir'-r^djg.p - ^7^^^7-/55^,^ ^ l_^>^'^^'^Pojg^^ I . (24) 

Putting Eqs. (l23l) and together, certain cancellations occur, and we find the 
final equation: 

r^r^dlSg^p = r'^r^dlSg^^ . (25) 
The latter has a solution 

5gaf3 = (1 + cu)hafi{x), (26) 

where c is an arbitrary constant to be fixed by the boundary condition^. The two 
sign choices considered above will hereafter be encoded in the value of c. We'll keep 
this constant unspecified till the end of our calculations. 
Using the solution fl26|) in equation fl2T|) we find: 

- ^ (□ V - V;.V./i^) + 2H^h^, - ^H%,h^ 

+H^c{h^,-%,h2) = T^u/2, (27) 

with its trace equation 

-3H'^ch = T/2. (28) 

Combining the above two equations, introducing the Lichnerowicz operator which 
in our case satisfies: 

AlH^, = -□ V + 8H^Ku - 2H%,hZ , (29) 

and using the standard techniques (see, [151 [IB] for recent discussions), we obtain 
the following expression for the perturbations: 

, ^ 1 ^ _1- 1 ^ 

Al-6H^ + 2H^{c+1) 3^^"-n-6i/2 + 2i72(c+i) 

O'c^'"' -□ - + 2m{c + 1) QH^c -□ - QH^ + 2H^{c + 1) ' ^ ^ 



^The expression in is not a most general solution of pS)) . however, it can be selected among 
all the solutions by specifying appropriate boundary conditions (see below). 
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The expression for the physical one-particle exchange amplitude reads as follows: 



" / '''^^ (5 ' ^) ^' -□ - 6H^ I 2HHc 

where T''^'^ denotes another conserved probe source. This should be compared with 
the amplitude for a massless graviton on dS space 

A = / d'xv^r^^h,^ = I d'^V^{T" \^^Q^J ,^ - l^' .a-QH^ ^} ' ^^^^ 

or with the amplitude for a massive graviton of mass M on dS space 

Am = j d'x^T"--h,,, = 

For c = — 1 the amplitude fl3Tl) is equivalent to that of a massless tensor field on dS 
in GR (l32l) . The solution for the background plus its perturbation in this case reads 
as follows: 

g^^ = {lT\u\)^^u + {lT\u\)h^u, (34) 

where h^^ is given in fl30l) . Note that these solutions corresponds to choosing ([8]) as 
the second boundary condition. 

For c > 2 one gets a massive graviton on the dS background [I?] and a massive 
scalar with the graviton mass = 2if^(c+ 1) and the scalar mass = 2H^{c + 
1) — GH"^. Moreover, the scalar couples to the stress-tensor with the 1/c suppressed 
strength as compared with the gravitational coupling. The metric for the solutions 
takes the form: 



gf,u = {It \u\)-ff,u + {^±c\u\)hf,^. (35) 

This solution corresponds to choosing in addition to ([T]) the following boundary 
condition: dug^,u\u=o+ = =F ((1 + c)%i. - cg^^). 

For c = the solution exist only for conformal sources with T = 0, for which 
one gets a special massive tensor on dS background with enhanced symmetry |18j . 

The boundary conditions with values of c other than c = — 1, c = and c > 2 
give rise to instabilities: for c < — 1 one gets a tachyonic tensor field (which implies 
that its helicity-0 component is a ghost) and a ghost-like scalar; for — 1 < c < one 
gets a massive tensor and a tachyonic scalar ghost; For < c < 2 one gets massive 
tensor and a tachyonic scalar. 
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Note that the last term in the expression for the field (!30|) is singular in the 
if = m — >■ limit. This term does not enter the linearized amplitude, but as it is 
well known, such terms typically give rise to strongly coupled behavior of massive 
theories [TH 1201 1211 [221 123]. This is a welcome feature in a classical theory as it 
provides a way to overcome the vDVZ discontinuity [21], as was first argued by 
Vainshtein p!9j (see also [20], and more recent works [25]). The magnitude of the 
strong scale should grow with c, as it's suggested by fISUl) . More detailed questions 
on its dependence on boundary terms and conditions are left open. The perturbative 
results obtained above have a limited applicabihty as the theory is expected to be 
strongly coupled. Moreover, perturbative stability is not a guarantee of a stability 
of the full nonlinear theory, however, it is a first and important step on the way to 
establish whether or not the theory could be viable. 



4. Hamiltonian 

In this section we derive the Hamiltonian for the theory ([2]). For this we use the 
standard ADM decomposition [B]: 

~goo = -{N'-Na''N^), ~9'' = jf2^ ^'' = ~9'' + ^^- (36) 

After somewhat lengthy algebra the additional term in the Lagrangian ([2]) can be 
written a^: 

(^J. -k')=V^ [N{ql - q^) - ^^^^ - 2qd^N^ , (37) 

where all indexes are raised by 7'-'; qij = \du'~iij = q = Y-'qij, and = duN^ ■ 

The expression in (137]) does not contain any time derivatives. Therefore, the 
canonical momenta in the extended theory ([2]) are the same as in GR. The Hamil- 
tonian density can straightforwardly be calculated 

= ^{NR^ + N,R^)5{u) T (^N{ql + q^) - ^'l^^' + 2Ar5„g^ + S .(38) 

Here S denotes the surface terms for both, the possible spatial boundaries, as well 
as the boundaries in the "w-dimension" 

S = 2Vj(7"'/'Ar,7r^--'X^i) ± 2m^du {^qN) . (39) 



^Note that in this section 7 refers to the 3D metric, as defined in ([M)) . while 7 denotes, as 
before, the 4D de Sitter metric. 
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The first two terms in (!38l) (the ones that are multiplied by S{u)) are those of GR 
with = -i?(3) + T-^TT^- - |7r2), and W = -2Vfc(7"^/^7r*^^), with being the 
canonical momenta of GR (see, e.g. [B], [3]). 

Since the additional terms in the Lagrangian ([2]) have no time derivatives fl37j) . 
the primary constraints of GR are preserved; the conjugate momenta for the lapse 
N and shift Nj are zero, P/v = Pnj = 0. Hence, variation of the Hamiltonian 6H 
under the variations of 6N (such that 6N{x,u)\u=±i = 0, and vanishing variation 
at u = and x = boundary) and 6Nj (such that 6Nj{x,u)\u=±i = and vanishing 
variation at m = at x = boundary )□ leads respectively to the following relations: 

R%iu) T 2m'^d^q T m'^ (gj + + = , (40) 



0. (41) 



Substituting these into the expression for the Hamiltonian fl38l) . one finds that the 
"bulk" terms all cancel and what is left is just the boundary terms: 



H{t) = j d^x j duHuix) = ±Am^ j d^x^/^{g'^ hj + g^'koi)\^\ 



(42) 



where we used the relation y/y{Nq + {NjduN^ /2N)) = \/g{g^^kij + g'^^koi), and 
dropped the surface term that appears in the GR Hamiltonian (that is the first 
term in 0391) ). Note that this is a Hamiltonian that follows from the Lagrangian 
([2]). If one adds additional surface terms as in (ITU]) , those terms should simply be 
subtracted from fj42|) to get the right Hamiltonian. 

For illustration we calculate the energy for the selfaccelerated solution ([9]) with 
a = 1 — The result is positive: 



H{t) = 6m^ J d^x^. (43) 

For the selfaccelerated solution with the growing a{u) in ^ the calculation of energy 
gives the same result (H3l) only after inclusion of the boundary term given in (fTOj) . 

As we see, the positive semi-definiteness of the Hamiltonian P2l) . in which the 
constraints (or algebraically determined relations) were used, depends on the bound- 
ary conditions in the u-direction. However, making these boundary terms positive 
semi-definite does not in general guarantee absence of instabilities, since the latter 
can be "hidden" in the constraint equations. One example of this is GR with a 
minimally coupled scalar of a negative kinetic term. The GR constraints put the 
Hamiltonian of this system to be zero, however, there are instabilities in the theory 
already at the classical level. 



^If the boundary is present in the x-space the GH boundary term should also be introduced. 
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In our case, the above derived results can be used to deduce the following im- 
portant observation: In the m — > limit the Hamiltonian fl42p goes to zero. This is 
in contrast with the 1 / term in the Hamiltonian for 4D massive gravity found by 
Boulware and Deser in jSj. Moreover, the expression fH2l) has no singular behavior 
in the field fluctuations, that was found in [8] as a source of various instabilities in 
massive gravity. Hence, even though there is no complete proof of the absence of 
instabilities in the full non-linear theory, the absence of the Boulware-Deser singular 
term in the Hamiltonian is a promising step forward. 

From Eq. (l40l) we find two equations for m = and u > respectively 

^RX=o = ±2m'V79\t- , VlduQ = (gj + + ^^^) > (44) 

where the rhs of the last equation is positive semi-definite. Similarly, we obtain from 
( HTl) the following equations for -u = and -u > respectively: 

Let us count the degrees of freedom. The variables A^|u=o and Nj\u=o can be fixed 
by gauge transformations. The variables duN\u=o and duNj\u=o can also be fixed 
after choosing the boundary conditions, for instance as A^|„=-|-i and Nj\u=±i, and 
using the equations flH|) and fHSl) . After fixing the boundary conditions what is 
left undetermined is the 12 variables ■yij\u=o,'^ij\u=o- Hence, in general this theory 
described 6 degrees of freedom, as we found it already in linearized calculations 
on the selfaccelerated background. For appropriate choice of boundary conditions 
these could be a massive graviton plus an additional scalar, which have no ghosts 
or tachyons, as it was shown in Section 3. For some particular boundary conditions 
though, due to enhanced symmetries of the linearized perturbations, the number of 
linear degrees of freedom gets reduced. In this case some of the equations in (jH]) 
and (H5|l should appear as constraints in the linearized theory. 

5. Discussions 

The extension of GR considered in this work is a convenient way of putting various 
theories of massive gravity in a single framework. All these theories, known in the 
linearized level, emerge as a consequence of choosing different boundary conditions 
in the auxiliary dimension. Moreover, the present framework provides a non-linear 
completion to these theories with the Hamiltonian that does not suffer from the 
problems found in Ref. [8]. The auxiliary dimension is just a convenient technical 
tool; it can in principle be "integrated out" entirely, and this should lead to GR 
amended by new terms in 4D. 

Most importantly, the extended theory admits the selfaccelerated solution with 
the spectrum of linear perturbations that has no ghosts or tachyons. In a general 
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case one obtains massive graviton and a scalar. This may have some cosmological 
signatures along the lines of Refs. [2S1EZ1- The vDVZ discontinuity of the linearized 
theory has to be overcome through the strong dynamics via the Vainshtein mech- 
anism [19] (see also [20]). If this is the case, then the theory is likely to have also 
short distance signatures [28], [29] . 

We end this section by a few comments. 

The auxiliary dimension discussed so far had a finite extent in the w-direction. 
It is straightforward to present a Lagrangian in which the w-direction is infinite: 

/ + 00 
duV^ {k% -k'-3) . (46) 
-oo 

The equations of motion of this Lagrangian have a selfaccelerated solution g^,^ = 
a{u)'j{x), where a{u) = e~", and as before, 7 denotes the 4D de Sitter metric with 
curvature R = 12m?. 

The Lagrangians ([2]) and ( H6l) . can be obtained by a certain truncation of a 5D 
theory. The 5D theory giving fH6l) can be defined as follows: 

^R + m,J dy^) {R{g) - R^{g^'^) - Sm^) j^^^^^^^^^^o , (47) 

where R^ is the 5D Ricci scalar, g'fj^ = {g^^,, 5(^5, 5(55}, A, 5 = 0, 1, 2, 3, 5, ?/ = u/rhc, 
m? = rricmc, and the substitutions in the last term are taken before the equations 
of motion are obtained, i.e., there is no variation w.r.t. (755 and g^^. To get the 
analogous expression for ([2|) one would have to drop the last term in the parenthesis, 
and set the integration w.r.t. y from — l/rric to +l/mc. The expression (1471) is 
somewhat similar to the DGP Lagrangian [5|, or its sign-flipped counterpart [TT] . 
with two crucial differences: (1) There is a subtraction of the R term from the R^ 
term in the bulk action; (2) There are no {55} or {/u5} equation^. 

Similar constructions with an auxiliary dimension can be considered for a scalar 
or vector, by adding the term — J du[{du(f))'^ + 0^ + ...], to the conventional scalar 
field Lagrangian, or the term —m^ j du[{duA^y' + ...] to the Maxwell Lagrangian 
(with a finite or an infinite range of integration). 
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^One gets back the analog of the {/i5} equation m Eq. © due to the Bianchi identities. It is 
really the absence of the {55} equation that makes a crucial difference. 
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